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SUMMARY 


A composite test specimen in the shape of a semicircular curved bar subjected to bending offers an excellent 
stress field for studying the open-mode delamination behavior of laminated composite materials. This is because the 
open-mode delamination nucleates at the midspan of the the curved bar. The classical anisotropic elasticity theory 
was used to construct a “multilayer” theory for the calculations of the stress and deformation fields induced in the 
multilayered composite semicircular curved bar subjected to end forces and end moments. The radial location and 
intensity of the open-mode delamination stress were calculated and were compared with the results obtained from 
the anisotropic continuum theory and from the finite element method. The multilayer theory gave more accurate 
predictions of the location and the intensity of the open-mode delamination stress than those calculated from the 
anisotropic continuum theory. 


NOMENCLATURE 


A,B,C,D 
C\ D 1 
A, B, C, D 


arbitrary constants associated with F for loading case of end forces P 

arbitrary constants associated with F for loading case of end moments M 

arbitrary constants associated with F for loading case of end moments M for isotropic materials 

inner radius of semicircular curved bar 

outer radius of ith layer of semicircular curved bar 

mean radius of semicircular curved bar, j(a + b ) 

outer radius of semicircular curved bar 

modulus of elasticity of single ply in fiber direction 

modulus of elasticity of single-ply transverse to fiber direction 

modulus of elasticity in r direction 

modulus of elasticity in 0 direction 

quadrilateral membrane element 

loading axis offset 

Airy stress function 

shear modulus of single ply 

shear modulus associated with r-6 system 

width of semicircular curved bar 

index associated with ith layer, i = 1,2,3 .... N 

anisotropic parameter, vf 
applied end moment 
total number of laminated layers 
applied end force 
radial distance 
radial location of gd 
radial location of ( a r )max 
radial location of (a^)max 



r o 

U r 

UQ 

P 

7 rO 

8 

Cr 

eo 

9 

VLT 

VtO, VOr 
Vtz.Vzt 
VQz,VzO 

&D 

G r 

( C^r)max 
( ^ r )max 
&8 
T t 9 
[ ] (0 
[ It 


radial location of zero gq 
displacement in r direction 
displacement in 9 direction 
rectangular Cartesian coordinates 

anisotropic parameter, ^J\ + 1 — 2 v t q) + 

shear strain in r-9 plane 

composite ply thickness 

strain in r direction 

strain in 9 direction 

tangential coordinate 

Poisson ratio of single-ply composite 

Poisson ratios 


delamination stress in C-coupon 
radial stress 

delamination stress for the case of end forces P, G r (r m , j) 
delamination stress for the case of end moments M> G r {r l m ) 
tangential stress 
shear stress 

quantity associated with ith layer 
quantity associated with ith layer 


1 INTRODUCTION 


One of the major causes of stiffness and strength degradations in laminated composite structures is the delam- 
inations between composite layers. In most engineering applications, laminated composite structures have certain 
curvatures (for example, curved panels and curved beams). If the curved composite structure is subjected to bending 
that tends to flatten the composite structure, tensile stresses can be generated in the thickness direction of the com- 
posites. Also, shear stresses could be induced if the bending is not a ''pure” bending. Under normal operations, if 
the above type of bending occurs cyclically, open-mode delaminations or shear-mode delaminations could nucleate 
at the sites of peak interlaminar tensile stresses or at the sites of peak interlaminar shear stresses. Continuation of 
these bending cyclings will cause the delamination zones to grow in size and ultimately cause the composite struc- 
tures to lose their structural integrity (loss of stiffness and strength) due to excessive delaminations. The type of 
delamination failure (open mode or shear mode) depends on which type of interlaminar strength (tensile or shear) is 
reached first. 

One of the most appealing geometries of a fatigue test coupon for studying the composite delamination phenom- 
enon is the semicircular curved bar shape (C-coupon). When such a test specimen is subjected to end forces (that is, 
nonpure bending), the peak radial stress (tension, if the bending tends to increase the radius of the curvature of the 
curved bar) and the peak shear stress induced in the curved bar will be identical in magnitudes but are out of phase 
in the tangential direction by £ Namely, the peak radial stress is located at the midspan point of the semicircular 
curved bar, but the peak shear stresses occur at both ends of the semicircular curved bar. The radial distance of both 
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the peak radial and the peak shear stresses are exactly the same. 1 The above nature of the semicircular curved bar of- 
fers an excellent situation for studying the initiation and subsequent propagation of delamination zones (open-mode 
or shear-mode) under cyclic loadings and for studying the fatigue behavior (degradation of stiffness and strength) of 
multilayered composite materials. In reference 1, Ko represented the multilayered composite semicircular bar with 
an equivalent continuous anisotropic material, and calculated radial locations and intensities of peak radial stresses 
induced in the curved bar subjected to end forces and end moments. Radial location and intensity of peak radial 
stress were calculated for different curved bar geometries (ratios of outer and inner radii) and for different degrees 
of anisotropy. Tolf 2 also conducted stress analysis of curved laminated beams using both continuous and discrete 
theories. He considered only the “pure” bending case. In this paper the multilayer theory (discrete theory instead of 
continuous theory) and the finite element method were used to perform similar delamination analysis of the multi- 
layered semicircular composite curved bar subjected to end forces and end moments. The resulting predictions of 
locations and intensities of peak radial stresses are compared with the results of the anisotropic continuum theory 
presented in reference 1 . 

2 COMPOSITE CURVED BAR 

Figure 1 shows the geometry of the composite curved bar (C -coupon) for delamination fatigue tests of composite 
materials. Because finite areas are needed for the load application points, both ends of the curved bar must be 
extended slightly. Thus, the C-coupon consists of a semicircular curved region with straight regions at both ends. 
Under the application of end forces P, the loading axis will have certain offset e from the vertical diameter of the 
semicircle. Thus, the loading condition on the C-coupon is the summation of the following two loading conditions 
(see fig. 2): (1) end forces P at the ends of the semicircle and (2) end moment M = Pe at the ends of the semicircle. 

Because the interface between 0° and 90° composite plies has the highest Poisson’s ratio mismatch in laying 
up the composite plies for fabricating the C-coupon, it is desirable to place the 90° or angle plies at the peak radial 
stress point to ensure that the delamination will nucleate at that point. Because of this demand, the precise location 
of the peak radial stress point must be known. The following sections will show how to determine the intensities 
and radial locations of peak radial stresses in the semicircular composite curved bar. 


3 ANISOTROPIC CONTINUUM THEORY 


For bending a linearly elastic continuous curved bar with cylindrical anisotropy, the Airy stress function P, 
written in cylindrical coordinate system, takes on the following functional forms: 3 

(a) for end forces P (fig. 2b) 

F = [At u< * + Br 1 ~^ + Cr + Dr In r] sin 6 (1) 

(b) for end moments M (fig. 2c) 

F = A'+ B'r 2 + C'r l+k + D'r X ~ k ( 2 ) 

where {A,B,C,D} and {A’ ,B' ,C' ,D'} are arbitrary constants that must be determined from boundary conditions, 
and the two anisotropic parameters (3 and k are respectively defined as 

(3) 

(4) 

For the isotropic case, /? = 2 and k = 1 . 
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The functional form given in equation (1) is also applicable to the isotropic case (/? = 2 ). However, the functional 
format given in equation (2) cannot be applied directly to the isotropic case by simply setting k = 1 . For the isotropic 
case, equation (2) must be expanded in the neighborhood of k = 1 using the relationship In 0 sa y4:(fc— l) _ ] 
where k - 1 < 1 to the following familiar form 4 

F = A + Br 2 + Cr 2 In r + D In r 1 51 


where {A,B,C,D} arc a different set of arbitrary constants. 

Stresses in the cylindrically anisotropic body may be expressed in terms of the stress function F as 

1 dF 1 d 2 F 



r r dr r 2 89 2 

(6) 


d 2 F 

dr 2 

(7) 


d 2 /F\ 
r ° 8r89 \ r ) 

(8) 

and the stress-strain relationships for the plane stress case are given by 



„ 1 VOr 

~ ~ Eo q ° 

(9)* 


VtB 1 

£> = - Ei' + 

(10)* 


1 

Ire = -p, — TtO 

Cr r$ 

(ID 

with the reciprocity relationship of 

VtB _ VBr 
E r ~ Eg 

02) 

The strains are related to the displacements through the following formulae: 



_ du T 

€t ~ ~d7 

(13) 


1 dug U r 

e 8 = — FF + — 

r dO r 

(14) 


1 / 1 du T dug ug \ 

tB 2 \r 89 + dr r ) 

(15) 


3.1 Stresses 


Substitution of equations (1) and (2) into equations (6) through (8) yields the following stress equations in terms 
ot the unknown arbitrary constants {A,B,D} or { B',C',D ' }: 

♦For the plane strain case, E r , Eg, u» r , and u r6 are simply replaced with — Mi , Ml gfc ( , , v r! v s ,\ 

v r „ /, UaV \ 1 -V„V„' 1 -V„U„ + Vf, ) md 

1 - v„v rt + ~u^)' respectively. 
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(a) for end forces P 


Ori. r, 0) = 

cre(r, 9) = 
T r fl( r,6) = - 


APr?-' - Bpr-P-' + — 


sin 9 


A0( 1 + P)r p ~' -Bp{ P)r P 1 + — 


Apr P- 1 - Bpr-P-' + — 


sin 6 


cos 6 


(b) for end moments M 


a r (r) = 2£' + C'(l + k)r k ~ x + D'{ 1 - k)r~ k ~ x 
og(r) = 2B' + C'k(l + AOr*' 1 - L>'k(l - k)r~ k ~ X 

T r8 = 0 


(16) 

(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 


Notice that the magnitudes of ay (equation (16)) and r r0 (equation (18)) for the case of end forces P are identical, 
but are out of phase in the 6 direction by ? . 


3.2 Displacements 


Using equations (9) through (12), (16) through (18), and (19) through (21), the displacements u r (equation (13)) 
and uq (equation (14)) may be integrated to give the following forms, neglecting the rigid body motion terms: 

(a) for end forces P 

uAr,S) = {Ar* [^ - O + /*>§*] + Br-> [± - (1 - 




ug(r. 


(b) for end moments M 


+ D( In r) 
Ar * 


-L-ttUsin!) 

E r Eg Ji 






+ D 


(In r) 


1 

E r E 6 


ii 

(i-t)]} cos< ’ 


u ,<r) = B'{2 r (^-^)} + C'{ ( l + «r‘(l^-f)} 

udr,S) - B '{ 2r “ ^)} 9 


3.3 Delamination Stresses and Their Locations 


( 22 ) 


(23) 


(24) 

(25) 


For the continuous (or single layer) curved bar, the two sets of unknown constants {A,B,D} and {B',C',D'} can 
be determined explicitly from the boundary conditions to give closed-form expressions for the stresses (equations 
(16) to (21)) and the displacements (equations (22) to (25)). 1,3 By using the extreme condition -^roy - 0, the 
functional expressions for the delamination stress (maximum value of a r ) and its radial location were derived in 
reference 1 for both of the aforementioned loading cases. 
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4 MULTILAYER THEORY 


Figure 3 shows the multilayer (AT -layers), semicircular curved bar subjected to both end forces P and end mo- 
ment M. The stress field and displacement field in each layer i(i = 1,2,... AT) for each loading case may be 
obtained from the results given in section 3. 

4.1 Boundary Conditions 

At each interface between layers i and * + 1(* = 1 , 2 , . . . N - 1) . the following boundary conditions for the 
continuities of stresses and displacements must hold (no sliding between layers). 


(a) for end forces P 



at r = <n: 

4'\ ai ,d) = ai i+l \a u d) 

(26) 


T$(a it d) = r r ( ’ +,) (a tI 0) 

(27) 


u^(a,,0) = o < i+1 >( o< ,0) 

(28) 


^(ai.fl) = u<‘ +1) (a t -,0) 

(29) 

The boundary conditions at the traction-free 
are 

inner surface (*—1 = 0) and outer surface ( * ' = 

N ) of the curved bar 

at r = oo = a: 

a^Ca.fl) = 0 

(30) 


T$g\a,6) = 0 

(31) 

at r = ax = b: 

4 N) (b,0) = 0 

(32) 


T^\b,6) = 0 

(33) 

(b) for end forces M 



at r = a*: 

a^Ca,-) = a^ +l) (a,) 

(34) 


Ur°U») = u^ +1) (o,) 

(35) 


«S°(0i) = 4 ,+ 1 >(a,) 

(36) 

at r = ao = o: 

a< 1} (o) = 0 

(37) 

at r = o n = b: 

aW(6) = 0 

(38) 
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As mentioned earlier, oy and t t q (equations (16) and (18)) for the case of end forces P have identical r depen- 
dency. Thus, if oy satisfies the boundary conditions, r T g will also satisfy the boundary conditions automatically. 
Therefore, the boundary conditions associated with t t q (equations (27), (31), and (33)) are not needed. 

For each loading case, each set of the previous boundary conditions will give 2+ 3( N - 1) = 3 jV - 1 equations 
for determining the 3 N unknowns A,y B,, Z?i( * = 1,2,... N ) for the case of end forces P, or 3 N unknowns 
B'i, C'i, D\{ i= 1,2,... TV) for the case of end moments M. 

The last equation needed for each loading case is the condition that the end force P or the end moment M is 
balanced by the stresses in the curved bar: 

(a) for end forces P 

-P = £f 4;\r,0)dr ; 0 = 0 (39) 

,=i • /o -i 

(b) for end moments M 

- M = [ ( r - r o )og(r)dr (40) 

i=i Ja '- 1 

where the negative signs in front of P and M are to increase the radius of curvature of the curved bar, and r 0 is the 
unknown radial location where ag = 0 . For pure bending we have 


/ cr g(r)dr = 0 

i=l •'*-» 


Therefore (since ro term vanishes), equation (40) becomes: 

N 


" r 0 * 

M = 52 / rag(r)dr 
i=i Ja '- 1 


(41) 


(42) 


4.2 Boundary Conditions in Final Forms 

After substitution of stress and displacement expressions given respectively in sections 3.1 and 3.2 into the 
boundary conditions given in section 4.1, the following final forms of the boundary conditions are obtained: 

(a) for end forces P 

for (equation (30)): 

A\P\aP' - + Di = 0 (43) 

for a (equation (26)): 

AiPia? - + Di - A l+ 1 1 af +1 + B i+ 1 ft + 1 a" A+1 - D i+ 1 = 0 ( 44) 

for (equation (32)): 


AnPn &N 1 — + Dm - 0 


(45) 
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for «<•> (equation (28)): 




,(*) 


(i + A)-dSrh 


E 


0i) (s* ~ 


+23, (In 


— 23,+ i(ln a,) 


( ' 

-d 

kJB (* +1) 

..(i+D \ 


1=0 

E ( g i+l) J 


.(i+l) 

Or 

n«+ 1) 


(46) 


for Ug'^ (equation (29)): 

w { ( 1 + A) - *4?} - W* A {( l - A) - *«} 


Eo 

+ D ‘^ (1 - * 4 r ^ 


-Ai+ 1 af ,+ 1 (i+ !) {( 1 + A+ 1 ) - ^flr + !> } + #«+ 1 a i A+1 {( 1 - A+ 1 ) - *4r + 1} } 

a 0 Eg 

(1 -»£**>) = 0 


(47) 


Equation (47) was obtained by taking the difference between the boundary conditions (28) and (29). This was done 
because the resulting expression (47) is simpler than using equation (29). 

for P (equation (39)): 


N 

E -“?-!> + £t(a~ A - a^>) + D,(ln a,- - In o,_i)l = P 

1=1 J 


(48) 


(b) for end moments M 
for cr^O (equation (37)): 


2B\ + C[( 1 + fcOo*-- 1 + DJ( 1 - fcOo-* 1 " 1 = 0 


for (equation (34)): 


(49) 


2 B' i+ c;(l+ kja?- 1 + D'( 1 - 
-2S; +1 - C' +1 (l + fc,+ i)a*‘* 1_1 - 23; +1 (l - Jfcf+Oor^- 1 = 0 


for <j( N ) (equation (38)): 


2B’ n + C' N ( 1 + k N )b k "~ l + D' n ( 1 - k N )b- kfl ~ 1 = 0 


(50) 


(51) 
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for u <‘> (equation (35)): 


4® 


+ CM(l + 


-d; jo -war*- (£]|o + i|))} 

B ‘*' { 2a ‘ “ 4 ^) } “ c ‘*‘ { (1 + fci * ,)o ‘"' (t^74 


_ ^gr 

(*+ 1) ^(l+l) 


+o< + i|a - 


—ki+l 


ki+ 1 ) Oi 


1 1 


U Jr 


ki + 1 fif 


(»+l) T ^(i+1) 


y g 

= 0 


(52) 


for ui 1 ^ (equation (36)): 


B'Alai 


E^ 


for M (equation (42)): 


- B ' i+ 1 


2oj 


e: 


(* + i) 1) 


= o 


(53) 


E [5,-(o? - o?_i) + c;fci(a * <+1 - o*iY ) - D'Ma -^ 1 - o^ 1 )] = — M 


(54) 


4.3 Delamination Stresses and Their Locations 


At exactly which layer the value of ay for each loading case will become maximum cannot be predicted until 
after all the unknown arbitrary constants (A,, B u C,} or {B[, C,-, D[} are determined from the appropriate boundary 
conditions given in section 4.2. Suppose (a r ) max (or ( o^Xm*), the maximum value of a r due to end forces P (or 
end moments M), occurs in the ith layer, then by using the extreme condition J^a r = 0, the radial location m (or 
r' m ) of ( 0y)max (or ( o^),™*) may be calculated from equation (16) (or equation (19)) as: 


(a) for end forces P 


(b) for end moments M 


— 


2 AifiAfii - 1) 


Pi. 

C< 


1 

TT 


(55) 


(56) 


And the delamination stresses ( a r ) max and ( a' ) max for the two loading cases may be written as 
(a) for end forces P 


(oy)max — ay 



(57) 
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(b) for end moments M 


(<4 W = a r (0 = 2B' i+ C'i(\ + kiHr'j 1 *- 1 + D'(l - k)(r' m )~* , ~ 1 


4.4 Delamination Stress in C-coupon 


(58) 


The delamination stresses ( a r )max (due to P) and ( a(.) max (due to M) do not occur at the same radial locations of 
the curved bar (see equations (55) and (56)). Thus, the delamination stress <jd in the C-coupon cannot be constructed 
by simply summing up ( oy)max and ( a' Xna* . The value of ctd must be evaluated at r = ru where the summation of 
the radial stress cr r ( r, f) due to P and the radial stress cr r ( r) due to M become maximum. Namely, 


Due to P Due to M 

= &r (r D , jJ + a r (r D ) 

where, the value of ru is calculated from the following extreme condition: 

Due to P Due to M 

^K r 'f) +a ' <r) ] =o 

which, after substitution of equations (16) and (19) and after performing differentiation, becomes 

AMfii ~ + B,A(A + l)r3 ft - Di 

+ (*, 2 - 1 )(C?rg + D'^) = 0 


(59) 


(60) 


(61) 


As will be seen later, the radial location of the delamination stress <jj) in the C-coupon is somewhere between 
r m and r' m (that is r m < r D < r' m ). 


5 NUMERICAL EXAMPLES 


The anisotropic continuum theory and the multilayer theory presented respectively in sections 3 (or ref. 1) and 
4 will now be applied to the delamination analysis of the composite C-coupon. One type of composite C-coupon 
under development has the following geometry and ply properties: 


inner radius 
outer radius 
loading axis offset 
width 

ply thickness 
mean radius 
radii ratio 


a = 2.1590 cm (0.85 in.) 
b = 2.9724 cm (1.17022 in.) 
e = 0.9525 cm (0.375 in.) 
h = 2.54 cm (1 in.) 

6 = 0.01506 cm (0.00593 in.) 
a m = ^ = 2.5657 cm (1.01011 in.) 

| = 1.3767 

a 

El = 17 .2369 x 10 10 N/m 2 (25 x 10 6 lb/in. 2 ) 
Et = 0.8274 x 10 10 N/m 2 ( 1.2 x 10 6 lb/in. 2 ) 
G lt = 0.4137 x 10 10 N/m 2 (0.6 x 10 6 lbAn. 2 ) 
vlt — 0.33 
vtl = 0.01584 
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The aforementioned C-coupon has 54 composite plies having the stacking sequence of [O 25 / + 15°/ — 15°/ — 
15°/+ 15°/0£]. 

5.1 Equivalent Continuum 

In order to apply the anisotropic continuum theory, the laminated composite C-coupon will 
be represented by an equivalent anisotropic continuum having the following effective material properties. 

Ee = 16.3220 x 10 10 N/m 2 (23.6731 x 10 6 lb/in. 2 ) 

E r = 0.8274 x 10 10 N/m 2 ( 1.2 x 10 6 lb/in 2 ) 

Gro = 0.4137 x 10 10 N/m 2 (0.6 x 10 6 lb/in. 2 ) 
v t q & 0 .01673 

Based on these effective material properties, the equivalent continuum representing the C-coupon has the values 
of anisotropic parameters as /3 = 7.7151 and k = 4.4416. 

5.2 Multilayer System 

For the purpose of applying the multilayer theory to the C-coupon, the extended linear regions at both ends will 
be neglected, and only the semicircular region subjected to two types of loadings (end forces P and end moments 
M , figure 2) will be considered. For simplification, each group of 25 layers of the 0° plies will be represented by one 
layer of anisotropic continuum, and the center region of 4 layers of ± 15 ° angle plies will be represented by another 
anisotropic continuum. Thus, the 54-layer composite will be represented by 3 layers of amsotropic continua. 

(a) 0 ° plies 

The inner ( i = 1) and the outer ( i = 3) layers have the following effective material properties: 

= E L = 17.2369 x 10 10 N/m 2 (25 x 10 6 lh/in. 2 ) 

£<•■> = Et = 0.8274 x 10 10 N/m 2 ( 1 .2 x 10 6 lb/in 2 ) 

G { $ = G lT = 0.4137 x 10 10 N/m 2 ( 0.6 x 10 6 lb/in. 2 ) 
v$ = 0 .01584 

which give the values of anisotropic parameters as pi = = 7.9272 and ki = ki = 4.5644. 

(b) ±15° plies 

The center layer (i = 2) has the following effective material properties: 

E ( g 2) = 4.8873 x 10 10 N/m 2 (7.0884 x 10 6 lb/in. 2 ) 

JS< 2 > = 0.8274 x 10 , 0 N/m 2 ( 1.2 x 10 6 lb/in. 2 ) 

G ( r f = 0.4137 x 10 10 N/m 2 (0.6 x 10 6 lb/in. 2 ) 
v ( T f ^ 0 .05590 

which give the values of the anisotropic parameters as = 4.2498 and kj = 2.4304. 

The 3 N (N = 3 ) boundary conditions for the aforementioned 3 -layer laminated system may be written in matrix 
forms (see appendix) for solving the 3 N unknown constants Ai, Bi, D t ( i = 1 , 2 , 3) , for the case of end forces P, 
or for solving the other set of 3 N unknown constants B' { , C - ,£>■(»' = 1,2,3), for the case of end moments M. 
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6 FINITE ELEMENT ANALYSIS 


To verify the solution accuracies of the analysis for which the actual loading condition of the 
C-coupon was represented with the superposition of two loading cases of the semicircular curved bar (see fig. 2), 
finite element stress analysis was performed on the semicircular curved bar (under two loading cases) and on the 
C-coupon using the structural performance and resizing (SPAR) finite element computer program. 5 Figures 4 and 5 
respectively show the SPAR finite element models set up for the semicircular curved bar and the C-coupon. Because 
of symmetry with respect to the x axis, only the half span of the semicircular curved bar and the C-coupon were 
modeled. Both systems were first reduced to three-layer systems as defined in section 5.2. Then layers 1, 2, and 
3 were respectively modeled in 10, 2, and 10 layers of quadrilateral membrane E4\ elements in the r direction. 
In the tangential direction, the quarter-circular region (0 < 0 < j)of the two systems was modeled with 90 E41 
elements. The extended region ( -e < x < 0) of the C-coupon was modeled with 25 layers of E41 elements in the 
x direction. 

The 6 = 90° plane for each model was allowed to move freely in the x direction (shown with rollers in figs. 4 
and 5) but not in the y direction. At the upper end of each model, only one point lying in the middle surface (og = 0 
point was found to be very close to the middle surface of the curved bar) was constrained to move freely in the y 
direction only (no movement in the x direction). Thus, the end of each model could rotate freely (shown with only 
one roller, figs. 4 and 5). The applied force P and the applied moment M were represented respectively with the 
distributions of T r g(r, 0)' and a$(r) obtained from the multilayer analysis. The sizes of the two SPAR models are 
listed below: 



Semicircular 



curved bar 

C-coupon 

JLOCs 

2093 

2668 

E4\ elements 

1980 

2530 


7 RESULTS 

7.1 Semicircular Curved Bar 


Figure 6 shows the distributions of a r in the 0 = f plane for the case of end forces P calculated from different 
theories. The values of (oy) ma x and r m calculated from different theories are indicated in the figure. The values of 
( °r)max calculated from different theories are quite close, except its location r m . The multilayer theory and the finite 
element method predicted close values of r m . The ( cr r ) max site predicted from the anisotropic continuum theory is 
located slightly closer to the middle surface than the ( oy) max sites predicted from the multilayer theory and SPAR. 
The ( o>)m a x site for the isotropic material is closest to the middle surface and is always located between the middle 
surface and the (oy)m a x site predicted from the anisotropic continuum theory. This can be seen more clearly from 
the plots shown in figure 7. 

The distance between the sites of ( cr r ) max predicted from the multilayer theory and the anisotropic continuum 
theory is 


( ^*fn) Anisotropic ( ^"m) Multilayer — (0 .4158 — 0 .3907)( b — fl) — 0 .0203cm (0 .0080in.) 

continuum 

which is 1.3554 times the single ply thickness of 0.01506 cm (0.00593 in.). 

Figure 8 shows the distributions of o r for the case of end moments M calculated from different theories. Unlike 
the previous case, the values of r' m and (cj') max calculated from different theories are quite close, showing that the 
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value of r f m is quite insensitive to the theory used. The multilayer theory predicted the shortest distance of r f m (that is, 
the ( CT^Xnax site is closest to the inner boundary of the curved bar). The site of ( o| )max predicted from the anisotropic 
continuum theory always lies between the middle surface and the (a(.)max site, based on isotropic theory. This is 
shown in figure 9. 

The distance between the sites of (a^) ma x predicted from the multilayer theory and the anisotropic continuum 
theory is given below 


( ^Anisotropic ” ( ^Multilayer = ( 0 .4346 - 0 .4327) ( b - a) = 0 .001545 cm ( 0 .000608 in.) 

continuum 

which is only 0.1025 times the single-ply thickness of 0.01506 cm (0.00593 in.), and is therefore, insignificant. 

Figures 10 and 1 1 respectively show the deformed shapes of the semicircular curved bar subjected to end forces 
P and end moments M. The dimensionless radial displacements ^ at 9 = 0 and 6 = j are shown in the figures. 

Table 1 summarizes all the values of (a r ) ma x. (oj-Xna*. r m , and r' m calculated from dif- 
ferent theories. 


Table 1. Intensities and locations of delamination stresses in semicircular curved bar. 


Theory 

End forces P 


End moments M 


Item 



Kb-a.1/ \ 

p V.c'rymax 

r m — Q 
b— a 

han 

u (°r)max 

r m~ Q 
h— a 

Anisotropic continuum 

1.4817 

0.4158 


1.5095 

0.4346 

Multilayer theory 

1.4864 

0.3907 


1.4988 

0.4327 

Isotropic continuum 

1.5135 

0.4205 


1.5222 

0.4338 

SPAR 

1.4816 

0.3935 


1.4979 

0.4398 


7.2 C-Coupon 

In order to determine the radial location rjj and magnitude of delamination stress a n for the C-coupon, the two 
radial stresses calculated from the semicircular curved bar due to P and M were summed up. The results are shown 
in figure 12 for multilayer, anisotropic, and isotropic cases. Notice that the radial location rj) predicted from all 
three theories lie between r m and r' m but closer to r m instead of r' m because the stress contribution due to P is larger 
than that due to M ( e < a m ) • The distance between the locations of au predicted from multilayer and anisotropic 
continuum theories is about 1.05 times the ply thickness. The finite element solution data points obtained from 
the C-coupon model lie in the vicinity of the two curves obtained from the multilayer and anisotropic continuum 
theories. The intensity and the radial location of <jd predicted from SPAR (C-coupon) are closest to those predicted 
from the multilayer theory. Table 2 summarizes the values of ctb and rjj predicted from different theories. 


Table 2. Intensities and location of 
delamination stresses in C-coupon 


Theory 

C-coupon under end forces P 

h(b—a) _ rn-o 

~P CT D 

Anisotropic continuum 

2.0415 

0.4212 

Multilayer theory 

2.0399 

0.4017 

Isotropic continuum 

2.0783 

0.4241 

SPAR (semicircular curved bar) 

2.0353 

0.3935 

SPAR (C-coupon) 

2.0405 

0.3935 
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Figure 13 shows the deformed shape of the C-coupon subjected to end forces P. The dimensionless radial 
displacement ^ at midspan and at the free end are shown in the figure. 

8 CONCLUSIONS 

The multilayer theory was developed for delamination analysis of a semicircular composite curved bar subjected 
to end forces and end moments. The difference between the radial locations of the delamination stress (maximum 
radial stress) predicted from the multilayer theory and from the anisotropic continuum theory was approximately 1.4 
times the ply thickness for the case of end forces and about 1/10 of the ply thickness for the case of end moments. 
The superposition method (namely, by summing up the two radial stresses induced in the semicircular curved bar 
subjected to end forces and end moments), used to construct the delamination stress in the C-coupon, gave reasonably 
accurate intensity of the delamination stress for the C-coupon. The finite element analysis of the C-coupon gave the 
radial location of the delamination stress in the C-coupon much closer to that predicted from the multilayer theory 
than from the anisotropic continuum theory. 


Ames Research Center 
Dryden Flight Research Facility 
National Aeronautics and Space Administration 
Edwards, California, December 14, 1988 
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Figure 4. Finite element model for multilayered semicircular curved bar under bending. 
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Figure 5. Finite element model for laminated composite delamination test C-coupon. 
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Figure 8. Distribution of radial stress, curved bar under end moments M\ b/a = 1.3767. 


M 



Figure 9. Plots of locations of (crj.) mai as a function of b/a for different values of k. 
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Figure 10. Deformed shape of curved bar under end force P. 
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Figure 12. Locations of dclamination stress an in C-coupon predicted from different theories 



Figure 13. Deformed shape of C-coupon under end forces P. 
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